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ABSTRACT 


A  flow  visualization  study  of  natural  convection  in  a 
liquid  filled  rectangular  enclosure  with  a  small  heater  pro 
truding  into  it  from  one  vertical  wall  was  conducted.  The 
top  and  bottom  horizontal  surfaces  of  the  enclosure  were 
heat  exchangers  whose  temperatures  could  be  varied  independ 
ently.  The  fluid  in  the  enclosure  was  water.  The  Baker 
electrochemical  technique  which  utilizes  a  pH  indicator  was 
used  for  flow  visualization.  Photographs  were  taken  of  the 
two  dimensional  flow  for  three  different  locations  of  the 
heater  on  the  vertical  wall.  Nusselt  numbers  for  each 
heater  location  and  for  a  range  of  Rayleigh  numbers  were 
also  determined.  Visual  results  indicate  the  presence  of 
a  dual-cell  configuration  within  the  enclosure.  There  is 
a  buoyancy  driven  upper  cell  and  a  shear  driven  lower  cell. 
The  flow  of  the  upper  cell  follows  the  geometry  of  the 
enclosure.  Though  the  uncertainty  inherent  in  the  experi¬ 
ment  was  high  data  calculations  suggest  a  trend  that  as  the 
heater  is  raised  within  the  enclosure  the  Nusselt  number 
decreases . 
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NOMENCLATURE 


Symbol 

Description 

Units 

A 

Area 

.2  „  2 
m  or  ft 

g 

Acceleration  of  gravity 

Ft/sec2 

Gr 

Grashof  number 

Dimensionless 

h 

Convection  heat  transfer 
coefficient 

Btu/hr-ft2-°F 

k 

Thermal  conductivity 

Btu/hr-ft-°F 

kf 

Fluid  thermal  conductivity 

Btu/hr-ft-°F 

kg 

Thermal  conductivity  of  plexiglas 

Btu/hr-ft-°F 

ki 

Thermal  conductivity  of  foam 
rubber  insulation 

Btu/hr-ft-°F 

kp 

Thermal  conductivity  of  phenolic 

Btu/hr-ft-°F 

kss 

Thermal  conductivity  of  304 
stainless  steel 

Btu/hr-ft-°F 

L 

Characteristic  length 

Inches 

Nu 

Nusselt  number 

Dimensionless 

Pr 

Prandtl  number 

Dimensionless 

^cond 

Heat  loss  via  conduction  through 
vertical  walls 

Btu/hr 

^conv 

Energy  convected  into  fluid 
from  heater 

Btu/hr 

^in 

Energy  into  cartridge  heater 

Btu/hr 

R 

Resistance  of  precision  resistor 

Ohms 

Ra 

Rayleigh  number 

Dimensionless 
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Ra 

Equivalent  thermal  resistance  to 
conduction  through  viewing  wall 

Hr-°F/Btu 

Rb 

Equivalent  thermal  resistance  to 
conduction  through  wall  opposite 
viewing  wall 

Hr-°F/Btu 

Rc 

Equivalent  thermal  resistance  to 
conduction  through  wall  holding 
heater 

Hr-°F/Btu 

Rcond 

Total  thermal  resistance  to 
conduction  with  respect  to  heat 
loss  through  vertical  walls 

Hr-°F/Btu 

Rj 

Individual  thermal  resistance 
to  conduction  (j=l,7) 

Hr-°F/Btu 

Ta 

Ambient  temperature 

°F 

Tf 

Average  temperature  of  heat 
exchangers 

°F 

Th 

Cartridge  heater  temperature 

°F 

Ti 

Thermocouple  temperatures  (i=l,8) 

°F 

Ts 

Average  temperature  of  heater 

°F 

V1 

Voltage  across  heater 

Volts 

V2 

Voltage  across  precision  resistor 

Volts 

W 

Uncertainty  of  variables 

Various 

Thermal  diffusivity 

2  , 

Ft  /sec 

s 

Coefficient  of  expansion 

1/°F 

V 

Kinematic  viscosity 

Ft2/sec 
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I .  BACKGROUND 


A.  INTRODUCTION:  STATEMENT  OF  PROBLEM 

The  rapid  advancements  in  the  field  of  miniaturization 
of  electronic  components  is  well-known.  Technology  has 
brought  us  from  the  large-scale  integration  (LSI)  device  to 
the  very  large-scale  integration  (VLSI)  device  in  just  a  few 
short  years.  Since  the  first  1  kilobit  RAM  semiconductor 
device  was  introduced  in  1970  the  number  of  memory  cells  on 
a  single  device  has  grown  to  64K  while  the  size  of  the  device 
has  decreased.  One  such  device,  the  64K  RAM  Micron  chip  pro¬ 
duced  by  Tiny  Micron  Technology,  Inc.  is  only  22,000  square 
mils,  or  approximately  the  size  of  one  of  these  printed 
letters  [1].  Yet  microminiaturization  advancements  are  far 
from  over  --  256K  and  greater  RAM  chips  are  on  the  immediate 
horizon. 

As  in  most  rapidly  developing  technologies  there  are  var¬ 
ious  problems  associated  with  the  advancements  that  have  been 
made  in  the  miniaturization  of  electronic  devices.  One  of 
the  primary  problems  that  has  evolved  is  that  of  suitably 
cooling  components  so  that  they  will  perform  reliably  over 
long  periods  of  time.  For  though  the  size  of  the  individual 
component  has  shrunk  considerably  its  junction  temperature 
has  remained  fairly  constant.  It  follows  that  the  decrease 
in  device  size  has  been  kin  to  an  increase  in  heat  flux. 


To  get  an  idea  of  the  magnitude  of  this  decreased  size/ 

increased  heat  flux  relationship  one  only  has  to  heed  the 

words  of  Chu  [2].  He  asserts  that  the  heat  flux  of  a  micro- 

circuit  chip  is  but  two  orders  of  magnitude  less  than  the 

S  7 

heat  flux  of  the  sun  (10  vice  10  Watts  per  meter  squared) . 
Further,  the  microcircuit  chip  must  dissipate  its  energy  at 
strict  maximum  junction  temperatures  of  between  125  to  150 
degrees  C.  The  integration  of  a  multitude  of  devices  on  a 
single  printed  circuit  board  leads  to  an  increase  in  the 
temperatures  in  which  the  device  operates.  The  high  heat 
flux  of  the  microcircuit  chip  coupled  with  the  increase  in 
environmental  temperature  very  often  leads  to  the  problem 
of  component  failure  resulting  from  inadequate  cooling.  It 
is  this  problem  that  has  recently  prompted  extensive  research 
into  the  techniques  of  cooling  microelectronic  devices. 

B.  DESIGN  CONSIDERATIONS:  PROPOSED  SOLUTIONS 

Many  papers  have  been  published  outlining  alternative 
methods  of  thermally  controlling  microelectronic  devices. 
Hannemann  [3]  describes  the  thermal  considerations  that  must 
be  taken  into  account  at  each  level  (subassembly,  assembly, 
and  system  levels)  of  the  device  design  process.  Kraus,  et 
al  [4]  assert  the  need  for  thermal  analysis  even  after 
design.  They  propose  a  strategy  that  yields  conservative 
estimates  of  the  heat  flux  associated  with  electronic  cir¬ 
cuitry.  Further,  Bar-Cohen,  et  al  [5]  indicate  that  numerical 
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analysis  of  the  heat  transfer  problem  should  be  relied  upon 
to  cope  with  the  analytically  intractable  nature  of  even 
standard  packaging  configurations. 

Once  the  initial  design  step  has  been  completed  a  reli¬ 
able  cooling  scheme  must  be  employed.  A  number  of  such 
schemes  including  conduction,  forced  air  cooling,  and  direct 
liquid  cooling  have  been  proposed  by  many  authors.  Bergles , 
et  al  [6]  contend  that  of  the  many  techniques  available, 
direct  liquid  cooling  may  be  the  most  desirable  mode  of  heat 
transfer  when  considerations  such  as  weight,  size,  or  cost 
are  taken  equally  into  account  with  the  heat  transfer  coef¬ 
ficient.  Simons  and  Moran  [7]  support  this  view  by  stating 
that  liquid  coolants  offer  the  capability  to  achieve  heat 
transfer  coefficients  one  to  two  orders  of  magnitude  higher 
than  that  of  air. 

Many  authors  have  offered  alternative  methods  of  utiliz¬ 
ing  immersion  cooling.  One  such  scheme  is  proposed  by  Kraus 
and  Bar-Cohen  [8]  who  analyze  submerged  condensers.  They 
describe  a  model  whereby  the  heat  from  the  immersed  elec¬ 
tronic  devices  is  removed  via  bubble  pumped  augmented  con¬ 
vection  interacting  at  the  submerged  condenser  surface. 

This  study  focuses  on  direct  liquid  cooling  by  natural 
convection.  More  specifically,  since  most  microelectronic 
devices  are  seated  on  printed  circuit  boards  packaged  very 
near  each  other,  investigation  is  limited  to  natural  convec¬ 
tion  in  a  small  liquid  filled  rectangular  enclosure.  One 
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advantage  of  natural  convection  cooling  is  that  there  are  no 
pumps,  pipes,  or  any  other  equipment  within  the  enclosure; 
the  manufacture  of  a  relatively  small  module  may  be  assured. 

C.  IMMERSION  COOLING:  ANALYTICAL  STUDIES 

Numerous  analytical  studies  of  natural  convection  in  rec 
tangular  enclosures  have  been  made  since  1966  when  Wilkes 
and  Churchill  [9]  initiated  finite  difference  analysis  for 
such  configurations.  This  work  was  followed  by  many  others 
including  Szekely  and  Todd  [10]  who,  in  1970,  presented  a 
finite  difference  analysis  of  natural  convection  in  a  rec¬ 
tangular  cavity  in  which  two  opposing  vertical  walls  provide 
the  temperature  gradient.  Comparing  their  numerical  results 
with  experiment  they  concluded  that  within  the  Rayleigh  num¬ 
ber  range  of  sixteen  hundred  to  sixty- thousand  the  thermal 
response  of  a  system  is  approximately  independent  of  the 
Rayleigh  number.  They  further  noted  a  trend  for  secondary 
cells  to  develop  once  the  Rayleigh  number  exceeded  three 
hundred  thousand. 

In  1975  Chu,  et  al  [11]  studied,  via  finite  difference 
methods  and  experiment,  the  natural  convection  in  a  rectan¬ 
gular  enclosure  in  which  the  heat  source  was  a  long,  hori¬ 
zontal  strip  located  n  one  vertical  wall.  Focusing  on 
heater  location,  1  ize,  and  enclosure  size  a  conclu¬ 

sion  was  reached  t  .  a  3  Rayleigh  number  increases  the 
heater  location  for  me.,  m -im  circulation  shifts  downward  but 
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that  heater  size  seems  to  have  little  effect  on  circulation. 


They  contend  that  secondary  cells  develop  at  the  upper  sur¬ 
face  as  the  height  of  the  enclosure  is  increased. 

Markatos  and  Pericleous  [12]  have  recently  published  a 
study  of  buoyancy-driven  flow  in  a  square  cavity  with  dif¬ 
ferentially  heated  sides.  Secondary  cells  were  also  observed 
in  their  study. 

D.  IMMERSION  COOLING:  EXPERIMENTAL  STUDIES 

Megerlin  and  Vingerhoet  [13]  collected  natural  convection 
data  using  fluorocarbon  dielectric  fluids  in  an  enclosure 
modeled  after  a  microelectronic  module.  One  significant  con¬ 
clusion  reached  was  that  cold  plate  orientation  has  a  signi¬ 
ficant  effect  on  the  heat  transfer  --a  fifty  percent  increase 
is  attained  using  a  horizontal  vice  vertical  cold  plate.  An 
experimental  investigation  using  an  undivided  and  partially 
divided  water  filled  two-dimensional  rectangular  enclosure 
was  performed  by  Nansteel  and  Greif  [14], 

In  1972  Baker  [15]  published  the  results  of  his  experi¬ 
mental  work  in  which  he  used  resistors  immersed  in  Freon. 

He  concluded  that  the  heat  dissipation  obtained  with  Freon 
is  at  least  four  times  greater  than  that  obtained  using  the 
same  assembly  with  stagnant  air.  Lin  and  Akins  [16]  investi¬ 
gated  natural  convection  heat  transfer  in  a  fluid  filled  cube 
where  all  six  wall  temperatures  were  subjected  to  a  step 
change  in  temperature.  Their  conclusion  was  that  a  variety  of 
circulation  patterns  evolve  depending  on  the  size  of  the  cube. 
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In  1980  Ostrach,  et  at  [17]  published  the  results  of  an 
experiment  that  used  small  particles  introduced  in  silicone 
oils  that  filled  various  sized  rectangular  enclosures  with 
height  to  length  ratios  less  than  unity.  Secondary  cells 
were  evident  in  their  study.  Finally,  an  experiment  was  per¬ 
formed  by  Bohn  and  Kirkpatrick  [18]  which  studied  high  Rayleigh 
number  natural  convection  in  a  cublical  enclosure  filled  with 
deionized  water  using  various  combinations  of  heated  walls. 

They  concluded  that  heating  from  below  along  with  simultaneous 
vertical  wall  heating  and  cooling  enhances  heat  transfer. 

The  above  discussion  is  not  intended  to  be  a  definitive 
synopsis  of  all  works  in  the  design  considerations  of  micro¬ 
electronic  devices,  nor  in  the  numerical  analysis  and  experi¬ 
ments  conducted  on  natural  convection  in  enclosures.  References 
cited  are  merely  intended  to  be  representative  of  the  works 
that  have  been  done  in  these  areas.  They  serve  to  provide  an 
excellent  point  of  departure  for  the  remainder  of  this  study. 

For  a  more  extensive  review  of  the  works  done  concerning  natu¬ 
ral  convection  in  rectangular  enclosures  the  reader  is  referred 
to  Catton  [19]. 

E.  OBJECTIVES 

The  objectives  of  this  study  are: 

1.  to  design  and  build  an  apparatus  that  approximates  a 
liquid  filled  microelectronic  package  with  appropriate  printed 
circuit  board  spacing  that  is  cooled  by  heat  exchangers  acting 
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as  the  top  and  bottom  horizontal  surfaces.  The  heat  source 
is  to  be  a  rectangular  protrusion  modeled  geometrically  after 
a  40-pin  dual-in-line  package  (DIP)  that  is  attached  to  and 
horizontally  traverses  the  length  of  one  vertical  wall. 

2.  to  observe  the  steady  state  natural  convection  flow 
in  the  apparatus  using  a  variety  of  heater  and  heat  sink 
temperatures  as  well  as  various  locations  of  the  heating 
element . 

3.  to  analyze  data  obtained  and  determine  Nusselt  and 
Rayleigh  numbers  for  each  apparatus  configuration. 

This  study  is  intended  to  serve  as  a  basis  for  future 
heat  transfer  experiments  in  natural  convection  using  vari¬ 
ous  enclosure  sizes  and  multiple  heated  protrusions. 
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II.  EXPERIMENT 


A.  GENERAL  DESIGN  CONSIDERATIONS 

The  flow  visualization  apparatus  (Fig.  1)  consists  of  a 
small  liquid  filled  rectangular  enclosure  with  a  heater 
modeled  geometrically  after  a  40-pin  DIP  protruding  into 
the  enclosure  from  one  of  the  vertical  walls.  This  con¬ 
figuration  is  meant  to  simulate  a  microelectronic  component 
on  a  printed  circuit  board.  The  top  and  bottom  horizontal 
surfaces  of  the  enclosure  are  such  that  their  temperatures 
can  be  varied  independently.  Heat  exchangers  (Fig.  2)  used 
in  conjunction  with  a  circulating  bath  act  as  these  horizon¬ 
tal  surface.  A  cartridge  heater  housed  in  a  stainless  steel 
cover  (Fig.  3)  acts  as  the  40-pin  DIP  model.  This  assembly 
is  hereafter  referred  to  as  the  heater.  A  calibrated  resis¬ 
tor  put  in  series  with  the  heater  and  its  power  supply  allows 
the  input  power  to  be  accurately  determined.  Eight  thermo¬ 
couples  are  placed  throughout  the  enclosure  for  temperature 
measurements.  Two  are  placed  on  each  heat  exchanger  and 
four  on  the  cartridge  heater  cover. 

Several  criteria  were  established  to  be  used  as  guidelines 
for  design.  The  following  is  a  list  of  these  criteria  and 
their  implications. 
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NOTC:  *U  OIWHSIOM  III  INCUTS  O  THIMWCOim  10CATIONS 


1 .  Dimens ions 

As  previously  noted,  the  heater  was  to  model  geomet¬ 
rically  a  40-pin  DIP.  Micrometer  measurements  of  the  height 
and  width  of  a  packaged  40-pin  DIP  were  made  and  the  rectan¬ 
gular  heater  cover  is  made  to  approximate  these  dimensions. 
Further,  the  enclosure's  height  was  to  be  6  times  the  verti¬ 
cal  dimension  of  the  heater  cover  and  its  width  was  to  be 
between  3  and  4  times  the  heater  cover's  horizontal  dimen¬ 
sion.  This  provides  a  good  approximation  to  the  spacing  of 
a  vertical  printed  circuit  board  placed  in  an  electronic 
module . 

2 .  Heater  Location 

The  apparatus  was  to  be  made  such  that  the  heater 
could  be  raised  or  lowered  on  the  vertical  wall  and  that  the 
width  of  the  enclosure  could  be  made  to  vary  from  1  to  2  inches. 
This  would  provide  a  wide  range  of  enclosure  sizes  and  heater 
locations  for  investigation.  It  was  determined  that  the 
ability  to  move  the  vertical  wall  on  which  the  heater  is 
attached  is  better  for  preserving  fluid-tight  integrity  and 
smoothness.  Machine  milled  interconnecting  parts  are  success¬ 
fully  used  to  meet  this  and  the  variable  chamber  width 
requirements . 

3.  Visualization  Technique 

The  thymol-blue  pH  indicator  flow  visualization  tech¬ 
nique  proposed  by  Baker  [20]  was  to  be  used.  This  is  a  non- 
intrusive  technique  that  uses  a  pH  indicator  in  distilled 
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water.  When  a  voltage  is  applied  across  the  fluid  it  percep¬ 
tibly  darkens  near  the  negative  electrode  due  to  a  local 
change  in  pH  induced  by  proton  transfer.  If  there  is  fluid 
motion  away  from  the  cathode,  as  there  is  in  this  experiment, 
the  darkened  fluid  traces  the  path  of  motion  before  dissipat¬ 
ing  back  to  the  original  lighter  color  and  neutral  pH  of  the 
bulk  fluid.  In  addition  to  its  non-intrusive  quality  no 
density  differences  occur  within  the  fluid  during  voltage 
application  since  the  pH  indicator  remains  an  ion  in  solu¬ 
tion.  The  specifics  of  the  fluid  preparation  are  discussed 
in  Chapter  III.  After  much  testing  it  was  determined  that 
the  use  of  this  technique  requires  any  metal  coming  into 
contact  with  the  liquid  solution  (e.g.,  heat  exchangers, 
heater)  to  be  relatively  inactive  with  respect  to  the 
Galvanic  series.  This  prevents  corrosion  as  indicated  by 
O'Conner  [21].  To  this  end,  a  stainless  steel  shim  strip 
running  the  length  of  the  vertical  wall  opposite  the  heater 
is  used  as  the  cathode.  The  heater  acts  as  the  anode. 

This  arrangement  provides  excellent  visualization  of  the 
two  dimensional  flow  of  the  darkened  solution. 

4 .  Other  Considerations 

One  eighth  inch  plexiglas  is  used  as  the  vertical 
viewing  wall  to  minimize  optical  distortion.  Further,  all 
vertical  walls  are  heavily  wrapped  with  foam  rubber  insula¬ 
tion  to  minimize  heat  loss  via  conduction  through  them 
during  experimental  operation. 
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B .  COMPONENTS 


1 .  Heater 

The  heater  assembly  (Fig.  3)  consists  of  an  8  inch 
long,  1/2  inch  wide,  1  inch  high,  and  3/16  inch  thick  304 
stainless  steel  channel  that  houses  a  narrow  rectangular 
cartridge  heater  backed  by  phenolic  insulation  that  runs 
the  length  of  the  channel.  The  cover  has  four  1/16  inch 
deep  slots  located  on  the  horizontal  faces  2.75  inches  from 
the  ends  of  the  channel.  In  each  of  these  slots  a  copper- 
constantan  thermocouple  is  spot  welded  nearest  the  cover's 
vertical  face.  The  thermocouple  wires  run  out  the  back  of 
the  slots  into  the  channel.  The  slots  are  filled  with  steel 
epoxy.  Three  evenly  spaced  screw  holes  are  tapped  into  the 
top  and  bottom  of  the  cover  so  that  the  entire  assembly  can 
be  attached  to  a  vertical  wall  in  the  enclosure.  Prior  to 
installation  the  cover  is  polished  to  a  mirror-like  finish 
using  180  to  600  grit  special  silicon  carbide  grinding  paper 
and  then  cleaned  with  common  rubbing  alcohol. 

2 .  Heat  Exchangers 

The  two  aluminum  heat  exchangers  measure  1.5  by  2.5 
by  11.0  inches.  The  top  heat  exchanger  is  seen  in  Fig.  2. 
They  each  have  three  rectangular  channels  separated  by  1/8 
inch  walls.  The  channels  are  machine  milled  into  a  1.375  by 
2.5  by  9.25  inches  solid  piece  of  stock  aluminum  and  then 
covered  with  a  1/8  inch  thick  plate  secured  by  ten  no.  4 
screws.  Slots  1/8  inches  deep  are  milled  8.0  inches  apart 


on  the  side  opposite  the  cover  plate  to  accept  vertical  walls. 
One  slot  is  3/16  inches  in  width  and  the  other  3/4  inches  in 
width.  Both  slots  traverse  the  width  of  the  heat  exchanger. 

In  addition,  two  1/16  inch  deep  slots  are  milled  2.75  inches 
from  the  outer  slots  and  they  run  to  within  1  inch  of  the 
other  side  of  the  exchanger.  These  last  two  slots  have 
copper-constantan  thermocouples  spot  welded  into  them  and 
are  covered  with  steel  epoxy.  A  1.5  by  8.0  inch  0.001  gage 
stainless  steel  shim  facing  is  epoxied  in  between  the  two 
outer  slots  on  the  side  of  the  exchanger  opposite  the  cover 
plate.  It  is  positioned  flush  against  the  corner  of  the 
heat  exchanger  that  is  one  inch  from  the  thermocouple  slots. 
This  facing  is  the  relatively  inactive  portion  of  the  ex¬ 
changer  that  acts  as  the  horizontal  wall  of  the  chamber. 

The  1.5  by  2.5  by  0.875  inch  headers  are  also  machined 
from  single  pieces  of  aluminum  stock.  They  each  have  a  1/4 
inch  od,  1/8  inch  id  exterior  pipe  for  coupling  with  a  cir¬ 
culating  bath  and  a  2  square  inch  hollow  interior  positioned 
so  as  to  get  maximum  flow  in  each  exchanger  channel.  The 
headers  are  held  to  the  main  body  of  the  exchangers  by  a 
single  no.  6  screw.  All  seams  are  fitted  with  epoxy  prior 
to  assembly. 

Both  heat  exchangers  were  tested  for  temperature 
distribution  and  leakage  Drior  to  use  in  the  experiment. 

They  were  coupled  to  a  circulating  bath  and  temperature 
measurements  were  taken  using  the  two  installed  thermocouples 
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and  an  additional  similar  thermocouple  that  was  randomly 
moved  about  the  exchanger  surface.  The  thermocouples  were 
run  to  a  digital  pyrometer  with  an  accuracy  of  1  degree 
Fahrenheit.  It  was  determined  from  this  test  that  the  heat 
exchangers  are,  upon  reaching  steady  state,  at  a  constant 
temperature  throughout;  i.e.,  the  pyrometer  reading  was  the 
same  for  each  thermocouple. 

3 .  Test  Chamber 

The  test  chamber  seen  in  Fig.  1  forms  a  rectangular 
enclosure  measuring  1.5  inches  in  width,  8.0  inches  in  depth 
and  6.0  inches  in  height.  The  top  and  bottom  horizontal  sur 
faces  of  the  enclosure  consist  of  the  stainless  steel  shim 
facing  of  the  heat  exchangers.  Three  of  the  vertical  walls 
are  entirely  plexiglas.  The  wall  holding  the  heater  is  10.5 
inches  wide,  24  inches  high,  and  11/16  inches  thick.  The 
heater  is  secured  to  the  center  of  this  wall  with  six  no.  4 
screws  and  rubber  o-rings  that  are  counter  sunk  into  the 
wall.  In  addition,  a  1/2  inch  diameter  hole  is  located  1 
inch  from  one  end  of  the  wall  to  accommodate  the  cartridge 
heater  power  leads  and  thermocouple  wires  from  the  heater 
cover.  The  wall  has  two  1/4  inch  deep  slots,  one  3/16 
inches  in  width  and  the  other  3/4  inches  in  width,  8  inches 
apart  and  14  inches  long  running  parallel  along  its  height. 
These  slots  receive  the  adjacent  vertical  walls.  They  also 
allow  the  wall  and  subsequently  the  heater  to  be  raised  or 
lowered  while  maintaining  fluid  tight  integrity. 
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The  other  two  plexiglas  vertical  walls  are  the  ends 
that  comprise  the  1.5  by  6.0  inch  dimensions  of  the  enclosure. 
They  are  both  2  inches  by  6.25  inches  but  one  is  11/16  inches 
thick  while  the  other  is  1/8  inches  thick.  The  thinner  wall 
is  used  for  visualization  and  photography  and  was  chosen  to 
be  thin  in  order  to  minimize  optical  distortion.  Both  walls 
slide  into  the  1/4  inch  slots  of  the  two  larger  adjacent  walls 
and  into  the  1/8  inch  deep  slots  on  the  heat  exchangers. 

The  final  vertical  wall  of  the  enclosure  is  a  10.0 
inch  by  6.0  inch  plexiglas  frame  that  house  a  glass  plate. 

The  plate  is  installed  with  epoxy  and  will  be  used  for  future 
studies  in  conjunction  with  laser  doppler  velocimetry.  Like 
the  other  large  wall  it  has  two  1/8  inch  deep  slots  that 
receive  the  adjacent  walls.  The  wall  has  lips  that  slide 
into  the  heat  exchanger  slots  as  well  as  two  1/16  inch  id 
standpipes  made  of  1/8  inch  stainless  steel  tubing.  One 
standpipe  is  located  near  an  upper  corner  and  the  other  is 
diagonally  opposite.  They  are  used  for  filling/emptying  the 
enclosure  and  relieving  internal  pressure  during  experimen¬ 
tal  operation. 

C.  ASSEMBLY 

Once  the  heater  has  been  attached  to  the  large  wall  the 
apparatus  is  assembled  with  the  exception  of  the  glass  wall 
framed  by  plexiglas.  Dow  Corning  732  RTV  adhesive/sealant 
is  applied  externally  to  all  seams  not  involved  with  visuali¬ 
zation,  to  the  power  lead/thermocouple  hole,  and  to  the  heads 
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of  all  heater  attaching  screws.  A  thin  coat  of  Dow  Corning 
3140  RTV  coating  and  silastic  is  used  on  all  internal  seams 
and  on  all  of  the  external  seams  of  the  smaller  walls.  The 
glass  wall  framed  by  plexiglas  is  then  fitted  into  place 
along  with  a  taut  1/4  inch  wide  strip  of  0.001  gage  302 
stainless  steel  shim.  The  strip  is  centered  internally  on 
the  wall  and  runs  its  height.  A  few  inches  of  the  strip 
traverse  the  seams  between  the  wall  and  the  heat  exchangers 
(both  top  and  bottom)  to  the  outside  of  the  test  chamber. 

This  strip  acts  as  the  cathode  for  the  thymol-blue  visuali¬ 
zation  technique.  The  sealants  are  then  applied  to  all 
remaining  seams  and  the  apparatus  is  tightened  by  specially 
made  clamps.  Fig.  4  is  a  photograph  of  the  assembled  apparatus. 

The  vertical  wall  of  the  apparatus  that  holds  the  heater 
has  to  be  removed,  placed  at  the  desired  heater  height,  and 
then  resealed  between  test  runs.  If  removing  the  wall  between 
test  runs  is  not  desired  the  apparatus  can  be  turned  upside- 
down  to  attain  a  different  heater  height. 

D.  INSTRUMENTATION 

1 .  Power  to  the  Heater 

The  heater  is  run  in  series  with  a  precision  resistor 
that  was  measured  to  have  a  2.0317  ohms  resistance  on  a 
Rosemount  commutating  galvonometer  bridge.  This  arrangement 
allows  a  voltage  measurement  across  the  resistor  to  be  made 
and,  using  its  known  resistance,  the  current  through  the 


29 


Figure  4.  Photograph  of  Experimental  Apparatus 
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circuit  can  be  calculated.  A  measurement  of  the  voltage 
across  the  heater  is  also  taken.  The  product  of  the  current 
and  the  voltage  across  the  heater  gives  input  power  to  the 
heater.  Both  voltages  are  measured  using  a  Dana  digital 
multimeter . 

Because  of  the  180  watt,  120  volt  rating  of  the  car¬ 
tridge  heater  two  different  power  supplies  are  used  in  order 
to  attain  certain  voltage  levels.  A  Lambda  regulated  DC 
power  supply  is  used  for  voltages  up  to  60  volts.  A  Variac 
voltage  AC  transformer  is  used  for  voltages  greater  than  60 
volts.  A  concern  about  using  the  latter  was  the  fluctuating 
voltage  provided  by  house  current  and  its  possible  effect  on 
measurements.  As  it  happened,  temperature  measurements  re¬ 
mained  stable.  HoweveT ,  voltage  measurements  across  the 
heater  fluctuated  up  to  3.0%.  Subsequently,  calculations 
of  input  power  to  the  heater  are  less  accurate  with  this 
mode  of  operation  than  the  power  calculations  made  when  the 
Lambda  regulator  is  used.  A  better  method  when  using  the 
Variac  would  have  been  to  put  it  in  series  with  a  line  volt¬ 
age  regulator  in  order  to  eliminate  voltage  fluctuations 
caused  by  house  current. 

2 .  Voltage  Across  the  Fluid 

The  voltage  source  for  the  electrochemical  technique 
is  also  a  Lambda  regulated  DC  power  supply.  Voltages  across 
the  electrodes  range  from  3.0  to  8.0  volts  depending  upon 
the  contrast  of  the  fluid  required  for  photography.  The 
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voltage  is  measured  with  the  Dana  digital  multimeter.  The 
stainless  steel  shim  strip  running  the  height  of  the  verti¬ 
cal  glass  wall  framed  by  plexiglas  acts  as  the  cathode  while 
the  stainless  steel  heater  cover  acts  as  the  anode.  Elec¬ 
trical  access  to  the  heater  cover  is  via  a  wire  spot  welded 
to  the  bottom  center  stainless  steel  attaching  screw  that 
holds  the  heater  assembly  to  the  vertical  wall  opposite  the 
cathode . 

Originally  the  heater  was  utilized  as  the  anode  and 
the  top  heat  exchanger  acted  as  the  cathode.  This  arrange¬ 
ment  proved  unacceptable  because  flow  paths  could  not  easily 
be  discerned  from  the  viewing  perspective.  The  inclusion  of 
the  shim  strip  as  the  cathode  solved  this  problem  by  limit¬ 
ing  the  darkened  flow  to  a  single  path  located  in  the  center 
of  the  enclosure. 

3.  Temperature  Measurement 

Eight,  30  gage  copper -const ant an  (teflon/teflon  in¬ 
sulation)  thermocouples  dispersed  throughout  the  enclosure 
are  used  for  temperature  measurements.  Two  thermocouples 
are  located  on  each  horizontal  face  of  the  heater  cover  and 
two  on  each  heat  exchanger.  All  thermocouples  were  fabricated 
using  a  Dynatech  welder  and  calibrated  on  a  Rosemount  commu¬ 
tating  bridge  in  conjunction  with  a  platinum  resistance  tem¬ 
perature  standard  and  calibration  bath. 

The  thermocouples  run  to  a  Newport  digital  pyrometer 
model  267  where  temperature  readings  are  made.  The  pyrometer 
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gives  readings  to  the  nearest  degree  Fahrenheit.  If,  during 
the  measuring  process,  a  reading  fluctuates  between  two 
temperature  extremes  it  is  assumed  to  be  the  average  of 
the  two.  Finally,  the  ambient  temperature  is  measured  with¬ 
in  two  feet  of  the  apparatus  with  another  thermocouple  run 


to  the  pyrometer. 


III.  EXPERIMENTAL  PROCEDURE 


A.  ELECTROCHEMICAL  FLUID  PREPARATION 

As  noted  earlier,  the  thymol -blue  pH  indicator  flow  vis¬ 
ualization  technique  as  described  by  Baker  [20]  is  used. 
Further,  O’Conner  suggests  using  thymol-blue  sodium  salt 
vice  pure  thymol-blue  (thymolsulphonephthalein)  crystals  in 
the  solution  to  avoid  filtering  of  undissolved  crystals  [21]. 
This  advice  is  followed  in  the  preparation  of  the  fluid. 
First,  0.2  grams  of  thymol-blue  sodium  salt  are  dissolved 
in  1.0  liter  of  distilled  water.  This  gives  a  0.021  by 
weight  light  colored  solution  with  no  undissolved  solids. 
Next,  0.1  normal  sodium  hydroxide  is  added  drop  by  drop  to 
the  solution  until  it  titrates  to  a  dark  blue  color.  The 
dark  color  is  an  indication  of  a  slightly  basic  pH.  Finally, 
the  solution  is  titrated  back  to  its  original  lighter  color 
by  adding  drop  by  drop  0.1  normal  hydrochloric  acid.  The 
fluid  is  now  slightly  acidic  and  very  close  to  its  neutral 
pH  value  of  7.  After  preparation,  a  small  amount  of  fluid 
is  poured  into  a  shallow  glass  container.  Wires  attached  to 
a  1.5  to  7.5  volt  dry  cell  battery  are  used  as  the  electrodes 
and  the  fluid  is  tested  for  local  color  change.  All  tests 
indicate  adequate  color  change. 
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B.  APPARATUS  PREPARATION 

After  the  sealants  are  cured,  1/4  inch  tygon  tubing  is 
attached  to  the  heat  exchangers  using  stainless  steel  hose 
clamps.  The  exchangers  are  then  interconnected  in  series  to 
a  refrigerated  circulating  bath  placed  six  feet  from  the 
apparatus.  The  top  heat  exchanger  is  placed  upstream  cf  the 
bottom  exchanger.  Common  tap  water  is  used  as  the  refriger¬ 
ated  circulating  coolant.  The  apparatus  is  then  laid  hori¬ 
zontally  on  its  largest  wall.  1.5  feet  of  1/8  inch  tygon 
tubing  is  attached  to  each  standpipe  and  the  enclosure  is 
filled  with  the  thymol-blue  solution  via  a  simple  siphon 
system.  Manual  manipulation  of  the  apparatus  and  standpipe 
tubing  allows  most  bubbles  to  be  driven  out  of  the  enclosure. 

A  specially  made  steel  leveling  stand  with  four  height 
adjustable  legs  and  a  slot  for  the  largest  vertical  wall  had 
been  made  to  hold  the  apparatus.  The  stand  is  placed  on  a 
piece  of  plywood  under  which  lies  two  inches  of  foam  rubber 
insulation.  This  system  serves  to  damp  out  any  vibrations 
occurring  in  the  laboratory.  The  apparatus  is  placed  verti¬ 
cally  upright  onto  the  leveling  stand  and  secured  in  place. 

The  two  larger  vertical  walls  are  then  layered  with  3/4  inch 
thick  insulation.  After  insuring  that  the  view  through  the 
test  chamber  to  a  white  background  is  unobstructed,  3/4  inch 
thick  insulation  is  placed  on  the  remaining  vertical  walls. 

The  apparatus  is  then  leveled  both  horizontally  and  vertically 
by  means  of  a  bubble  level  placed  on  the  top  heat  exchanger 
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in  two  orthogonal  directions  and  by  concurrent  adjustment  of 
the  stand  legs.  This  leveling  process  is  crucial  for  ensur¬ 
ing  horizontal  and  vertical  alignment  of  the  chamber.  Fig.  5 
is  a  photograph  of  the  assembled  apparatus  and  instrumentation. 

C.  TEST  PROCEDURE 

1 .  Heater  Location 

Test  runs  were  made  for  three  heater  locations:  (1) 
the  top  horizontal  face  of  the  heater  located  1.0  inch  from 
the  top  of  the  enclosure,  (2)  the  heater  centered  on  the 
vertical  wall  of  the  enclosure,  and  (3)  the  bottom  horizontal 
face  of  the  heater  located  1.0  inch  from  the  bottom  of  the 
enclosure . 

2 .  Initial  Instrument  Settings 

Temperature  measurements  of  all  thermocouples  are 
made  to  insure  that  the  apparatus  and  its  liquid  contents 
are  at  a  single  temperature.  If  this  criteria  is  met  the 
circulating  bath  to  the  heat  exchangers  is  turned  on  to  a 
predetermined  temperature.  Shortly  thereafter  the  power 
supply  to  the  heater  is  turned  on  and  the  voltage  slowly 
increased  to  another  predetermined  value  that  results  in 
the  approximate  desired  temperature  difference  between  the 
exchangers  and  the  heater.  The  values  of  the  settings  for 
the  circulating  bath  temperature  and  the  voltage  to  the 
heater  are  determined  soley  from  experience  garnered  from 
a  multitude  of  test  runs. 
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Instrument  Readings 


Thermocouple  temperature  measurements  and  voltage 
measurements  across  the  heater  and  resistor  are  taken  every 
twenty  minutes  until  successive  temperature  measurements 
remain  unchanged  to  within  the  accuracy  of  the  pyrometer. 
Measurements  are  then  taken  every  ten  minutes.  It  is  assumed 
that  steady  state  is  attained  within  the  enclosure  if  no 
thermocouple  measurement  changes  by  more  than  one  degree 
during  any  ten  minute  period.  This  process  takes  between 
two  to  four  hours  depending  on  the  initial  temperature  dif¬ 
ference  between  the  heater  and  the  exchangers.  The  smaller 
temperature  difference  takes  the  least  amount  of  time. 

Once  steady  state  is  attained  the  power  supply  that 
creates  the  voltage  across  the  electrochemical  fluid  is 
turned  on.  Using  a  digital  multimeter  as  an  indicator,  the 
voltage  is  increased  from  0.0  to  3.0  volts.  The  insulation 
is  then  removed  from  the  two  smaller  vertical  walls  and  the 
fluid  voltage  is  varied  until  a  visually  suitable  path  of 
darkened  fluid  ensues. 

4 .  Photographic  Techniques 

A  Nikon  Ftn  with  a  50  mm  f2  lens  is  used  for  all  black 
and  white  photography.  A  Canon  Ftb  with  a  50  mm  fl.4  lens 
is  used  for  all  color  photography.  Black  and  white  photo¬ 
graphs  are  taken  with  Kodak  Tri-x  pan  asa  400  film;  color 
photographs  are  taken  with  Kodak  Ektachrome  asa  64  daylight 
film.  Black  and  white  photographs  are  taken  first.  Both 


cameras  are  secured  to  a  tripod.  The  Nikon  is  placed  between 
14  and  15  inches  from  the  1/8  inch  thick  plexiglas  vertical 
wall.  A  combination  of  f-stop  settings  of  8  and  11  coupled 
with  shutter  speeds  of  1/60  and  1/125  of  a  second  are  used. 

The  Canon  is  placed  at  a  distance  of  between  12  and  13  inches 
from  the  viewing  wall.  Combinations  of  f-stop  2  and  4  with 
shutter  speeds  of  1/4  and  1/15  of  a  second  are  used  with 
the  Canon. 

Two  L-shaped  pieces  of  cardboard  are  attached  together 
so  as  to  create  a  rectangular  opening  approximately  the  size 
of  the  small  vertical  wall.  This  is  placed  between  the  camera 
and  the  enclosure.  A  cloth  is  hung  from  this  cardboard  aper¬ 
ture  and  draped  over  the  camera  during  photography  to  elimi¬ 
nate  camera  reflection  on  the  plexiglas.  Another  larger 
piece  of  cardboard  covered  with  white  paper  is  placed  1  foot 
from  the  rear  small  vertical  wall  to  act  as  a  backdrop. 

Finally  a  650  watt  colortran  multi-6  set  lamp  is  positioned 
2  feet  from  the  white  cardboard  backdrop  and  1.5  feet  from 
the  apparatus.  The  lamp  is  directed  towards  the  white  back¬ 
drop  and  turned  on  only  when  photographing  the  flow. 

A  major  concern  was  tnat  the  heat  generated  by  the 
lamp  would  disturb  the  flow.  However,  due  to  the  insulation 
remaining  between  the  lamp  and  the  enclosure  and  the  short 
interval  of  time  the  lamp  is  on  the  viewed  flow  remains  two 
dimensional.  This  scheme  provides  an  excellent  view  down  the 
length  of  the  enclosure  to  a  white  background  that  adequately 
contrasts  the  darkened  fluid. 
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D.  DATA  ANALYSIS 

There  are  six  length  dimensions  within  the  enclosure,  any 
one  of  which  could  be  chosen  as  the  characteristic  length, 

L,  for  determining  the  Nusselt  number  and  Rayleigh  numbers. 

The  1.0  inch  vertical  face  of  the  heater  is  chosen  since  it 
is  the  surface  that  not  only  adds  the  majority  of  energy  to 
the  system  but  it  also  is  oriented  in  the  direction  of  the 
gravitational  driving  force  of  free  convection.  The  proper¬ 
ties  of  the  fluid  are  evaluated  at  the  average  of  the  hot 
surface  (heater)  and  cold  surface  (heat  exchangers)  temperatures. 

1 .  Determination  of  Nusselt  Number 

The  Nusselt  number  is  defined  as  follows : 

Nu  =  (1) 


To  evaluate  the  heat  transfer  coefficient  an  estimate  has  to 
be  made  of  the  heat  lost  through  the  vertical  walls  of  the 
enclosure.  A  parallel  thermal  resistance  circuit  for  con¬ 
duction  through  the  three  vertical  walls  in  contact  with  the 
heater  is  employed  for  this  purpose.  Heat  loss  through  the 
vertical  wall  opposite  the  heater,  all  contact  resistances, 
and  conduction  via  th*1  heater  assembly  attaching  screws  are 
assumed  negligible.  Knowing  the  total  thermal  resistance 
for  conduction  the  heat  loss  is  calculated  using  a  variation 
of  Fourier's  Law  of  Conduction: 


Q 


cond 


I 


cond 


(2a) 
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where 

R  =  L 

Kcond 


(2b) 


The  length  of  the  conduction  path,  L,  for  all  materials  is 
known.  The  thermal  conductivity,  k,  for  each  material  is 
determined  using  a  table  of  properties  [22].  The  area.  A, 
is  assumed  to  be  that  of  the  cartridge  heater  perpendicular 
to  the  conduction  path.  Finally,  the  temperature  difference 
is  assumed  to  be  the  difference  between  the  ambient,  T_ , 

a 

and  the  heater  cover  temperature,  Tg ;  this  further  neglects 
the  convective  thermal  resistance  on  the  outside  of  the  in¬ 
sulation.  In  the  event  Tc  is  less  than  T  ,  heat  loss  is 

j  a 

assumed  to  be  zero. 

After  calculating  the  energy  lost  via  conduction  and 
having  already  determined  the  energy  into  the  system  (Qin) 
from  the  voltage  measurements  across  the  heater  and  the  pre¬ 
cision  resistor  a  simple  energy  balance  is  used  to  determine 
the  energy  convected  into  the  fluid  (Qconv) : 

^conv  *  ^in  "  ^cond 


where 


Q. 

xm 


V1  V2 


(3b) 


and 

Vj  *  voltage  across  heater 
V2  ■  voltage  across  resistor 
R  =*  2.0317  Ohms 
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Then,  using  Newton's  Law  of  Cooling, 


Q  =  h  A  (T 
xconv  v  s 


-  V 


(4) 


h  is  calculated.  T£  is  chosen  to  be  the  temperature  of  the 
heat  exchangers  which  are  the  ultimate  heat  sink  for  the 
enclosure.  The  exterior  surface  area  of  the  heater  cover 
is  A.  With  all  variables  in  equation  (1)  now  known  the 
Nusselt  number  is  calculated. 

To  ensure  the  accuracy  of  this  method  of  analysis  a  single 
iteration  is  performed  on  the  Nusselt  number  by  use  of  another 
thermal  resistance  for  conduction  circuit.  This  time  a  single 
circuit  with  a  lone  resistance  is  run  through  the  steel  heater 
cover  to  the  cartridge  heater  with  the  intent  of  determining 
the  cartridge  heater  temperature.  Using  the  heater  cover 
temperature  (Ts ) ,  length  of  the  conduction  path  through  the 
heater  cover  (L) ,  the  thermal  conductivity  for  304  stainless 
steel  (kss) ,  the  exterior  area  of  the  heater  cover  (A) ,  and 
the  previously  calculated  convection  energy  into  the  enclo¬ 
sure  (Qconv) »  it  is  a  simple  matter  to  determine  the  cartridge 
heater  temperature  (T^)  from  a  rearrangement  of  Fourier's  Law 
of  Conduction: 


T 


h 


conv 

nr~ 


+ 


(5) 
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Ts  in  equation  (2a)  is  then  set  equal  to  and  a  single 
iteration  is  performed  to  re-evaluate  the  Nusselt  number. 
Iterative  analysis  of  Nusselt  numbers  for  all  test  runs 
prove  to  fall  within  3.0%  of  the  initial  calculation  thus 
eliminating  the  necessity  for  further  iterations.  A  sample 
calculation  is  contained  in  Appendix  A. 

2 .  Determination  of  Rayleigh  Number 

One  definition  of  the  Rayleigh  number  is  as  follows: 


where 


Ra  =  (Gr)  (Pr) 


Gr 


Pr 


g  6  l  (T 


V 


v 

a 


(6a) 

(6b) 

(6c) 


In  the  calculation  of  the  Grashof  number  (Gr) ,  the  same 
characteristic  length  (L)  is  used  as  that  used  to  determine 
the  Nusselt  number;  that  is,  the  1.0  inch  vertical  face  of 
the  heater  cover.  In  addition,  the  same  temperature  differ¬ 
ence  is  used  as  that  in  equation  (4)  (i.e.,  the  average 
heater  cover  temperature  minus  the  average  heat  exchanger 
temperature) .  The  remaining  terms  in  equation  (6b)  are  eval 
uated  at  the  average  temperature  of  all  thermocouple  measure 
ments  by  a  cubic  polynomial  that  had  been  fit  to  data  from 
a  table  of  properties  [23]  for  temperatures  of  water  ranging 
from  50  to  100  degrees  Fahrenheit. 


The  Prandtl  number  (Pr)  is  also  evaluated  by  a  cubic 
polynomial  that  had  been  fit  to  data  from  the  same  table 
of  properties,  at  the  same  temperature,  and  for  the  same 
working  fluid.  The  Rayleigh  number  is  then  calculated. 

A  sample  calculation  is  contained  in  Appendix  A.  A  sample 
uncertainty  calculation  is  contained  in  Appendix  B. 


44 


IV.  RESULTS 


A.  VISUAL 

Photographs  of  the  natural  convection  flow  within  the 
enclosure  for  a  representative  number  of  heater  locations 
and  Rayleigh  numbers  are  presented  in  Figs.  6,  7,  and  8. 

Figure  6  depicts  flows  with  the  top  horizontal  face  of  the 
heater  located  1.0  inches  from  the  top  of  the  enclosure. 

Figure  7  depicts  flows  with  the  heater  centered  in  the  en¬ 
closure.  Finally,  Figure  8  depicts  flows  with  the  bottom 
horizontal  face  of  the  heater  located  1.0  inches  from  the 
bottom  of  the  enclosure. 

Some  general  observations  can  be  made  about  the  natural 
convection  flow  within  the  enclosure.  First,  the  flow  is 
dual-celled.  There  is  a  buoyancy  driven  upper  cell  and  a 
shear  driven  lower  cell.  The  location  of  the  cell  separa¬ 
tion  layer  is  dependent  upon  the  heater  location  within  the 
enclosure;  as  the  heater  is  raised  so,  too,  is  the  layer. 

Secondly,  the  flow  of  the  buoyant  driven  upper  cell  fol¬ 
lows  the  geometry  of  the  enclosure.  The  fluid  comes  off 
the  vertical  face  of  the  heater,  up  the  vertical  wall  to 
which  the  heater  is  attached,  across  the  cold  heat  exchangers, 
then  down  the  vertical  wall  opposite  the  heater.  Finally, 
the  fluid  moves  along  the  cell  separation  layer  where  it 
interacts  viscously  with  the  lower  cell  and  back  towards  the 
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6 


Figure  6.  Photographs  of  Flow  with  Heater  Ne 


heater.  No  definitive  comments  can  be  made  about  the  total 
motion  of  the  fluid  in  the  lower  cell;  only  a  very  slow  moving 
flow  was  observed  slightly  below  the  separation  layer.  This 
flow  was  crossing  the  width  of  the  enclosure  towards  the  ver¬ 
tical  wall  holding  the  heater. 

The  photographs  also  show  flow  being  generated  into  each 
cell.  This  is  the  result  of  a  recurring  problem  encountered 
using  the  Baker  flow  visualization  technique.  Baker  suggests 
applying  a  maximum  of  5.0  volts  across  the  fluid  when  elec¬ 
trode  spacing  is  0.5  cm.  These  limitations  preclude  genera¬ 
tion  of  bubbles  near  the  cathode.  Utilizing  the  limits  imposed 
and  a  network  of  cathodes  the  ensuing  flow  can  be  seen  as  it 
moves  away  from  the  area  of  interest.  In  this  experiment  a 
single  cathode  is  used.  The  area  of  interest  is  the  entire 
enclosure.  Hence,  it  is  necessary  to  apply  up  to  8.0  volts 
in  order  to  observe  the  flow  throughout  the  entire  enclosure. 
Bubbles  are  therefore  generated  that  cause  convective  mass 
transfer  up  the  stainless  steel  shim  strip.  The  buoyant  flow 
circulation  of  the  upper  cell  collides  with  the  mass  transfer 
associated  with  the  convected  bubbles  resulting  in  the  intru¬ 
sion  of  flow  into  the  cells.  Efforts  were  made  to  minimize 
this  effect  by  careful  adjustment  of  the  voltage  during  test 
runs . 

The  distinct  difference  between  the  photographs  of  Fig.  7 
clearly  illustrates  the  advantage  of  utilizing  a  stainless 
steel  shim  strip  as  the  cathode  for  the  visualization  technique 
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As  previously  noted,  the  heater  was  originally  utilized  as  the 
anode  while  a  heat  exchanger  acted  as  the  cathode.  This  pro¬ 
vided  an  excellent  darkening  of  the  fluid  but  also  proved  in¬ 
adequate  in  terras  of  flow  visualization  (right  photograph  of 
Fig.  7).  The  inclusion  of  the  shim  strip  as  the  cathode  allows 
visualization  of  the  two  dimensional  flow  that  occurs  in  the 
middle  of  the  enclosure  (left  photograph  of  Fig.  7). 

B.  QUANTITATIVE 

Table  I  summarizes  the  results  of  calculations  using  data 
obtained  for  each  test  run.  Based  on  these  calculations  the 
Nusselt  and  Rayleigh  numbers  were  plotted  on  a  log- log  graph 
(Fig.  9).  Since  the  temperature  measurement  uncertainty  in¬ 
herent  in  the  experiment  produces  Nusselt  and  Rayleigh  number 
uncertainties  ranging  from  five  to  twenty  percent  (see  Appendix 
II)  it  is  difficult  to  draw  definitive  conclusions  from  the 
data.  However,  there  does  appear  to  be  a  trend  suggesting 
that  as  the  heater  is  raised  within  the  enclosure  the  Nusselt 


number  decreases . 


TABLE  I 


Summary  of  Results 


Tg  (0F) 

T £  (°F) 

Qin  (Btu/hr ) 

Nu 

Ra  (x 

Heater 

Near 

Top 

65.75 

54.00 

26.41 

4.95 

0.99 

70. 2S 

50.00 

55.75 

6.12 

1.73 

72.25 

38.88 

102.68 

6.84 

2.17 

81.25 

38.88 

146.85 

7.65 

3.61 

82.50 

39.25 

145.52 

7.38 

3.85 

89. 7S 

36.75 

197.58* 

8.19 

5.33 

Heater 

Centered 

66. 7S 

57.00 

18.94 

4.23 

0.92 

69.25 

59.00 

25.99 

5.58 

1.07 

68. SO 

57.50 

28.79 

5.76 

1.09 

69.00 

50.50 

51.55 

6.12 

1.55 

66.75 

42.00 

68.15 

6.12 

1.49 

70.50 

40.63 

99.01 

7.38 

1.94 

70.50 

39.25 

103.80 

7.38 

1.95 

77.75 

39.75 

144.89 

8.46 

5.00 

75.50 

36.75 

147.37 

8.46 

2.61 

79.00 

39.50 

147.33 

8.28 

3.21 

80.00 

59.25 

148.06 

8.01 

3.38 

88.00 

40.00 

185.42* 

8.46 

5.01 

99.50 

38.00 

248.91* 

8.73 

7.95 

Heater 

Near 

Bottom 

68.25 

56.50 

28.78 

5.40 

1.13 

68.00 

51.00 

51.24 

6.66 

1.40 

68.50 

40.00 

98.12 

7.74 

1.70 

75.75 

40.00 

139.54 

8.64 

2.68 

Note : 

(1) 

*  indicates 

Variac  used  as 

power  supply 
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V.  CONCLUSIONS 


The  author  has  found  no  previous  natural  convection  study 
of  a  liquid  filled  rectangular  enclosure  with  a  heated  pro¬ 
trusion  attached  to  a  vertical  wall,  temperature  variant 
horizontal  walls,  and  the  remaining  walls  adiabatic.  However, 
some  investigators  using  similar  configurations  have  reached 
the  same  conclusions  about  to  be  presented. 

A.  VISUAL 

Two  cells  occur  within  the  enclosure  during  the  natural 
convection  flow.  The  upper  cell  is  driven  by  buoyant  forces 
and  essentially  follows  the  geometry  of  the  enclosure.  The 
lower  cell  is  shear  driven.  The  location  of  the  cell  separa¬ 
tion  layer  is  dependent  upon  the  height  of  the  heater  in  the 
enclosure . 

B.  QUANTITATIVE 

Calculations  indicate  a  trend  for  the  Nusselt  number  to 
decrease  as  the  height  of  the  heater  within  the  enclosure  is 
increased.  It  is  interesting  to  note  that  for  the  Rayleigh 
numbers  calculated  in  this  experiment  (on  the  order  of  105 
to  10 6 )  Chu,  et  al  [11]  plot  this  same  tendency  for  an  air 
filled  square  enclosure  that  has  a  heater  imbedded  in  a  ver¬ 
tical  wall  and  cooled  horizontal  walls.  It  is  Chu's  inves¬ 
tigation  that  most  closely  resembles  this  experiment. 
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VI.  RECOMMENDATIONS 


A.  IMPROVEMENT  TO  EXPERIMENT 

1 .  Apparatus 

Fluid  leakage  into  the  heater  assembly  occurred  during 
the  course  of  the  experiment.  This  caused  the  cartridge  heater 
to  fail.  Replacement  of  the  cartridge  heater  necessitated  dis¬ 
mantling  the  apparatus.  To  avoid  the  problem  of  cartridge 
heater  malfunction,  either  a  liquid  immersible  cartridge  heater 
should  be  used  or  else  an  improved  design  of  the  heater  cover 
should  be  sought. 

Though  the  photographs  proved  satisfactory  it  is  fur¬ 
ther  suggested  that  experimentation  with  various  camera  filters 
be  done.  The  contrast  of  the  observed  flow  may  improve  to  the 
point  where  excessive  voltage  across  the  fluid  may  not  be  re¬ 
quired,  thereby,  eliminating  the  bubble  generation  effect. 

2 .  Instrumentation 

In  order  to  lower  the  experimental  uncertainty  it  is 
suggested  that  a  more  accurate  temperature  measuring  instru¬ 
ment  be  used  for  thermocouple  readings.  Further,  to  preclude 
the  use  of  the  Variac  transformer  a  cartridge  heater  with  the 
same  power  (189  watts)  but  rated  at  60  instead  of  120  volts 
should  be  used.  This  will  allow  all  test  runs  to  be  made 
with  the  Lambda  power  supply  as  the  power  supply  to  the 
heater . 
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B.  ADDITIONAL  EXPERIMENTAL  WORK 

It  is  suggested  that  the  following  areas  of  study  be 
experimentally  explored: 

1.  Decrease  the  enclosure  width  to  see  if  any  additional 
cells  develop  within  the  enclosure. 

2.  Observe  the  effects  of  independently  varying  the  heat 
exchanger  temperatures.  In  conjunction  with  this,  do  not 
allow  the  coolant  temperature  to  go  below  fifty  degrees  Fah¬ 
renheit  to  avoid  the  density  inversion  that  occurs  in  water 
at  forty  degrees.  As  can  be  seen  in  Table  I,  the  heat  ex¬ 
changer  temperatures  for  some  test  runs  fell  below  this  value. 

3.  Continue  to  use  the  stainless  steel  shim  strip  as  the 
cathode  but  seek  methods  of  decreasing  the  effect  of  the 
bubbles  generated  when  voltage  is  applied  across  the  fluid. 
Alternatively,  design  a  network  of  cathodes  (as  suggested  by 
Baker)  that  encircles  the  enclosure. 
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APPENDIX  A:  SAMPLE  CALCULATIONS 

A.  DETERMINATION  OF  INPUT  POWER 

The  dimensions  of  the  cartridge  heater  are  8.0  by  0.25 
by  0.725  inches.  Table  II  is  the  data  that  was  recorded 
for  a  representative  run  on  3  November  1983  with  the  heater 
located  near  the  bottom  of  the  enclosure.  Refer  to  Fig.  1 
for  thermocouple  locations. 


TABLE  II 

Test  Run  Data  for  3  November  1983 

h  Ii  h  tjl  h  I±  h  li 

40  40  40  40  73  77  77  76  71  58.51  1.42 

Notes:  (1)  all  temperatures  in  degrees  Fahrenheit 

(2)  all  voltages  in  volts 


Using  equation  3b  the  input  power  is  calculated  to  be: 

Q.  *  (5-8  -  ^i-i1  (3.4123  Btu/hr-Watt)  =  139.54  Btu/hr 

in  (2.0317  ohms) 


B.  NUSSELT  NUMBER 

The  heat  loss  via  conduction  through  the  vertical  walls 
is  determined  by  constructing  a  thermal  resistance  network 
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(Fig.  10).  The  thermal  conductivity  of  each  material 
(Table  III)  is  determined  from  a  table  of  properties  [22]. 


RU8BER 

INSULATION  PLEXIGLAS 


Figure  10.  Thermal  Resistance  Network  for  Calculation  of  Losses 


TABLE  III 


Thermal  Conductivity  of  Materials 


Material 

K  (Btu/hr- 

rubber  insulation 

0.0225 

plexiglas 

0.0821 

phenolic 

0.0173 

304  steel 

8.8693 

Table  IV  is  a  summary  of  the  individual  resistances  for  Fig. 
10  that  were  determined  using  equation  2b. 


TABLE  IV 

Thermal  Resistances  to  Conduction 


Resistance 

Path (1) 

Material 

Value  (°F-h: 

R1 

a 

plexiglas 

116.9 

R2 

a 

rubber  insulation 

2562.7 

R3 

b 

plexiglas 

643.2 

R4 

b 

rubber  insulation 

2562.7 

R5 

c 

phenolic 

8.7 

R6 

c 

plexiglas 

20.1 

r7 

c 

rubber  insulation 

80.1 

Note:  (1)  conduction  path  through 

a.  viewing  wall 

b.  wall  opposite  viewing 

c.  wall  to  which  heater  is  attached 
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The  equivalent 


Therefore , 


resistance  for 

R  =  2679.6 
a 

Rb  =  3205.9 
Rc  =  108.8 


each  path  is: 
°F-hr/Btu 
°F-hr/Btu 
°F-hr/Btu 


Rcond 


tRaj  (Rb^  CRc> 
(Ra)(Rb)+(Ra)(Rc)+(Rb)(Rc) 


101.28  °F-hr/Btu 


Employing  equation  2a, 


Q 


cond 


(75,75  -  40.0) 
(101.28) 


0.047  Btu/hr 


The  Nusselt  number  is  determined  by  first  using  equation 
3a  to  calculate  the  energy  convected  into  the  fluid: 

Q  =  (139.54  -  0.047)  *  139.495  Btu/hr 
^conv 

Then,  the  heat  transfer  coefficient  is  calculated  as  (equa¬ 
tion  4) 

h  «  (139,495)  ^144)  -  35.10  Btu/hr-ft2 -°F 

(75.75-40.0) (16) 

The  thermal  conductivity  of  water  is  determined  from  a  table 
of  properties  [23] , 
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k£  =■  0.3383  Btu/hr-ft-°F 


Finally,  the  Nusselt  number  can  be  calculated  (equation  1) , 


Nu  =  i1/12)  =  8.64 

(0.3383) 


An  iterative  check  is  performed  by  determining  the  car¬ 
tridge  heater  temperature  (equation  5)  as 

T  =  IlZiilL-IllLiiil  +  75.75  =  77.96  °F 
n  (8.659)  (1.33) 


Reapplication  of  equations  2a,  3a,  4,  and  1  respectively  with 
Tg  set  equal  to  yields 


Q 


cond 


(77.9.6  -  71.0) 
(101.285) 


0.0687  Btu/hr 


QCOnv  =  (139.54  -  0.0687)  =  139.474  Btu/hr 

h  -  — (139.474)  (144) —  ,  35.10  Btu/hr-ft2-°F 
(77.96  -  71.0) (16) 


From  the  table  of  properties  [23]  the  fluid  conductivity 
remains 

kf  -  0.3383  Btu/hr/ft/°F 
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Therefore , 

Nu  .  (,35,:10j..(l/U)  ,  8  64 

(0.3383) 

The  percent  difference  in  Nusselt  numbers  after  one  itera¬ 
tion  is  zero  percent. 

C.  RAYLEIGH  NUMBER 

Using  the  property  table  [23]  at  the  average  temperature 
of  the  heater  and  heat  exchangers  the  Prandtl  number  and 
Grashof  number  can  be  calculated  as : 

Pr  =  8.322 

££  -  15.58  x  106  l/°F/cubic  ft 
v 

Gr  *  (15.58]  (77.96  -  40.0)  (1/1728)  =  322.4  x  103 

Therefore,  the  Rayleigh  number  can  be  determined  using 
equation  6a: 

Ra  -  (322.4  x  103)  (8.322)  *  2.683  x  106 
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APPENDIX  B:  UNCERTAINTY  ANALYSIS 


Uncertainties  in  the  experiment  are  evaluated  by  employing 
the  root  mean  square  method  of  Kline  and  McClintock  [24]. 
Table  V  is  a  summary  of  measurement  and  property  values  and 
their  associated  uncertainties: 


TABLE  V 


Uncertainty  in  Variables 


Variable 

gB/v2 

Uncertainty 

1/250 

kf 

1/300 

k 

5.0  % 

g 

k. 

1 

7.0  % 

k 

7.0  * 

P 

k 

5.0  % 

ss 

L  (1) 

0.001  in 

Pt 

1/600 

R 

0.0001  ohms 

T  (2) 

1.0  F 

V1(3) 

0.01  volts 

V,(4) 

2.0  volts 

v2 

0.01  volts 

Basis  for  Uncertainty 
assume  last  digit 
assume  last  digit 
[22] 

[22] 

[22] 

[22] 

resolution  of  measurement  device 
assume  last  digit 
calibration  on  galvonometer 
resolution  of  measurement  device 
resolution  of  measurement  device 
resolution  of  measurement  device 
resolution  of  measurement  device 


Notes:  (1)  all  dimension  measurements 

(2)  all  thermocouple  temperature  measurements 

(3)  Lambda  regulator  as  power  supply 

(4)  Variac  as  power  supply 
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The  following  variable  uncertainties  are  determined  using 
the  aforementioned  RMS  equation  on  the  data  presented  in 
Table  II  of  Appendix  A. 

2  2  2  ^/2 

(WV,)(V7)  (WV?)(V.)  (WR)(V.)(V2) 

WQ.  =  [( - i - £_)  +  ( - £ - ±_)  ♦  ( - i - i-)  ] 

in  R  R  R 

(0.01)  (1.42)j2.  ((0.01)  (58.51)}  2.  (  (0.0001)  (58.51)  (1.42)  } 
(2.0317)  (2.0317)  (2 . 0317) 2 


=  0.983  Btu/hr 


WQ 


cond 


WT  2 
[(—  )  +  ( 


WT  L 
—  )  +  (  — )  ♦  ( 


cond 


cond 


(W*cond>  <Ta 
(RconcP 


-  TJ 


1/2 


-)  ] 


[[(2)(l/2))2 

(101.28) 


{■(2)  (1/2)  j  +  ((0.0000S)(2))a] 
(101.28)  (101. 28) 2 


=  6.017  Btu/hr 


WQ 


conv 


tWcond’2  *  (W<5in)2l 

[ (0 . 017) 2  ♦  (0 . 983)  2  ) 


1/2 

1/2 


»  0.983  Btu/hr 
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1/2 


[((Wh)(L).  2  +  ((WL)(h)}2  +  f  CWkf )  (h)  (L)  ^  2  j 
(kf)  (kf)  (kf)2 


[  c  C2 . 95 )  (0.083))  2+  c(Q.001)(35.10)j  2  ^(0.001)  (35.10)  (0.083)^  ] 
(0.338)  (0.338)  (0.338)2 


0.727 


WRa  =  [((WGr)  (Pr))z  +  ((WPr)(Gr)K] 


=  [((21.0  x  103)(8.322))2  +  ((0.0) (322.0  x  105) )  ] 


=  178.85  x  10- 


The  above  uncertainty  calculations  aTe  intended  to  be 
representative  of  the  method  used;  calculations  performed  on 
every  variable  are  not  included.  A  summary  of  quantitative 
results  for  Nusselt  and  Rayleigh  numbers  and  their  respective 
uncertainties  in  relation  to  heater  temperature  minus  ex¬ 
changer  temperature  for  all  test  runs  is  given  in  Table  VI. 
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TABLE  VI 


Uncertainty  in  Experiment 


VTf 

Nu  $ 

Uncertainty  Ra 

(x  106) 

%  Uncertainty 

Heater 

Near  Top 

11.75 

4.95 

18.20 

0.99 

17.40 

20.25 

6.12 

11.73 

1.78 

10.52 

33 . 38 

6.84 

8.69 

2  .17 

7.00 

42.38 

7.65 

7  .85 

3.61 

5.95 

43.25 

7.38 

7.79 

3.85 

5.88 

53.00* 

8.19 

7 . 89 

5.33 

5.23 

Heater 

Centered 

9.75 

4.23 

21.57 

0.92 

20.83 

10.25 

5.58 

20.55 

1.07 

19.85 

11.00 

5.76 

19.28 

1.09 

18.53 

18.50 

6.12 

12.53 

1.55 

11.40 

24.75 

6.12 

10.26 

1 . 49 

8.86 

29.88 

7.58 

9.19 

1.94 

7.61 

31.25 

7.38 

8.97 

1.95 

7.35 

38.00 

8.46 

8.19 

3.00 

6.39 

38.75 

8.46 

8.13 

2.61 

6.31 

39.50 

8.28 

8.07 

3.21 

6.23 

40.75 

8.01 

7.97 

3.38 

6.10 

48.00* 

8.46 

8.07 

5.01 

5.52 

61.50* 

8.73 

7. S3 

7.95 

4.87 

Heater 

Near  Bottom 

11.75 

5.40 

18.19 

1.13 

17.40 

17.00 

6.66 

13.37 

1.40 

12.31 

28.50 

7.74 

9.43 

1.70 

7.90 

35.75 

8.64 

8.41 

2.68 

6.67 

Note : 

(1)  *  Indicates  Variac  used  as  power 

supply . 
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